Stability, causality and quasinormal modes of cosmic strings and cylinders 
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In this work we consider the evolution of a massive scalar field in cylindrically symmetric space- 
times. Quasinormal modes have been calculated for static and rotating cosmic cylinders. We found 
unstable modes in some cases. Rotating as well as static cosmic strings, i.e., without regular interior 
solutions, do not display quasinormal oscillation modes. We conclude that rotating cosmic cylinder 
space-times that present closed time-like curves are unstable against scalar perturbations. 
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(N ; I. INTRODUCTION 

Cylindrically symmetric space-times, cosmic strings and cosmic cylinders in particular, are an important subclass 
of the full axial geometries, with a symmetry structure in many ways more complex than the spherically symmetric 
$_i | backgrounds. Although not compatible with asymptotic flatness, cylindrical solutions have been studied in many 
different physically meaningful contexts. Definitions and structural aspects have been investigated in [H, 0. In 
cosmology and astrophysics they have been used to model the formation of structure , gravitational lensing effect 
i—i , [3] and gamma ray bursts Q. In condensed matter, they appear in the study of superconductivity and topological 
> 1 defects 0. 

Beside formalism issues and phenomenological applications, cosmic strings display a variety of different phenomena. 
They are important for the study of Einstein General Relativity in extreme situations, where the limits of the theory 
\ can be explored and information about very new phenomena can be gathered. For example, they have been used to 
. model time-machines 

@, S [13 ■ Cylindrical space-times frequently have cuts in the angular variable. In case there 
g ; is a rotation, one can imagine a voyager along the angular direction with an effectively superluminal velocity in a 
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closed space-time path, hence a time machine. Thus, such geometries are generally associated with the possibility of 
traveling back in time, always a puzzling possibility [Til . H2i . Il3j . Hence the question whether these space-times are 



physically feasible. Effective superluminal velocities are not necessarily a pathology. For instance, in [14j, it has been 
shown that in braneworlds a gravitational sign can be sent through the bulk in a way that it can reach a distant point 
in the brane faster than light propagating in the brane. However, such signs seem to be irrelevant, because most of 
5_j them advance, compared to a light sign, by negligibly small amounts. 

If we take causality as a fundamental principle, it is reasonable to conjecture that time machines cannot be con- 
structed in the real world, although they are possible as exact solutions of Einstein equations. It is expected that these 
time machines do not survive as physical solutions, but we do not yet know any physical mechanism that prohibits 
the existence of such backgrounds. 

In this work, we propose such mechanism, associating closed time-like curves to instabilities in certain cylindrical 
geometries. Taking the spherically symmetric cases as reference, not much has been stated about the stability of 
cosmic strings. This is specially true for rotating cosmic strings, which are relevant for time machine theoretical 
proposals. One exception, in the braneworld context, is the well known Gregory-Laflamme [l5j instability in the black 
strings. 

We have calculated the evolution of a massive scalar field propagating in different cosmic string and cosmic cylinder 
backgrounds. The relevant boundary conditions to be satisfied by a massive scalar field are presented. These conditions 
will be part of the quasinormal mode definition, which will be calculated considering these backgrounds. The static 
cases will be treated in the first place, because the rotating geometries present closed time-like curves (CTCs) and in 
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these cases additional analysis is necessary to render the problem tractable. Analyzing the behavior of the field and 
their quasinormal modes we could present evidence of the stability of some of these space-times. 

II. COSMIC STRING AND CYLINDER SPACE-TIMES 

The metric of a stationary and cylindrically symmetric space-time is given by an expression of the form 

ds 2 = —F dt 2 + 2M dtdcf) + L d<j) 2 + H dr 2 + S dz 2 , (1) 



where (F, M, L, H, S) are functions only of the radial coordinate r. All space-times studied here are described by this 
Ansatz. 

We start with the static cosmic cylinder background. It was obtained by Gott and Hiscock [J, [f| and its exterior 
solution is given by 

ds 2 = -dt 2 + dr 2 + dz 2 + (1 - 4/i) r 2 d<t> 2 , (2) 

where the coordinates range are — oo < t < oo, fo sin(#M) (1 — 4/z) 1 < r < oo, < <f> < 2ir, and — oo < z < oo. The 
interior solution is given by 

ds 2 = -dt 2 + dr 2 + dz 2 + f 2 sin(r / f ) 2 d(j> 2 , (3) 

with — oo < t < oo, < r < Tq6m, < <p < 27r, — oo < z < oo and fi is the linear mass density. The matching point 
is defined by r = r\, in the exterior radial coordinate and the corresponding interior radial coordinate is r = r s . They 
are related by 

f sin(6» M ) _ Q ... 

n = — , r s = r o 0M ■ (4) 

1 — 4/1 



The other constants are related as 

1 

4 L ~ ' " ' r 8vrf, 



t 1 = ~ [1 - COs(e M )} , P = ZZ=2 ' ( 5 ) 



where p is the energy density. The static cosmic string space-time can be obtained taking the limit fo — ► 0. For 
both static cosmic string and cosmic cylinder, we have assumed /i < 1/4. Additional details are presented in Q. The 
restriction on /x leads to 

= ■ (6) 

r 2 

The rotating cosmic cylinder solution was obtained by Jensen and Soleng (Io| . The exterior solution reads 

ds 2 = -dt 2 + dr 2 + dz 2 - %Jd<pdt + [(1 -4/i) 2 (Y + r ) 2 - 16 J 2 ] d(j) 2 , (7) 

where \x is the linear mass density and J the rate of angular momentum per unit length. The constant rg determines 
the origin of the exterior radial coordinate. For 

this space-time exhibits closed time-like curves. 

Inside the cylinder there are at least two possible regular solutions. The first scenario considered in this work is the 
"flower pot" model, 



ds 2 



2 



16J' . ., 
1 + — 7-(ri-r*) 
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8Jr 

dt 2 —d(j)dt + r 2 d(j) 2 +dr 2 + dz 2 , (9) 



where r s is related to the exterior parameter ro by 
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Here there is no closed time-like curve. 

The second type of internal solution treated is the so-called "ballpoint pen" model. It describes gravity generated 
by the following distribution of matter, 



Snp = A + n 2 



= a\(r s - r) , 



(11) 



where a < 1 is an arbitrary constant which can be adjusted such that we may or may not have closed time-like curves. 
This interior solution is described by 



ds 2 = -dt 2 - 2Md<t>dt + sm( ^ r ) 2 _ M 2^2 + dr 2 + dz 2 



(12) 



M = 2a 



i \ i f\ \ sin (VAr) 

(r - r s ) cos(V Xrj -= h r s 

v A 



(13) 



Imposing that exterior and interior solutions match continuously on the boundary r s of the rotating cosmic cylinder, 
the parameters are related by 



± 



VI- (1-4/^ 
(1 — 4^) arccos (1 — 4^) 



VA = 



arccos (1 — 4/i) 



(14) 



1 — cos(v / Ar s ) 



J = 



2a 



sin 



(VAr s 



(15) 



III. EVOLUTION OF A MASSIVE SCALAR FIELD 



The evolution of a minimally coupled massive scalar field f in a curved space-time is given by the equation 

□*(t,r,0,z) = fl 2 y(t,r,e,z) , (16) 
where (3 is the mass of the field and the differential operator □ has the form 

1 
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and g is the determinant of the metric. In the cylindrically symmetric space-time given by metric ([T]) we have 

1 d r [X HS ) / a ji, \ f Q x n\ 

- Ld 2 ^ + Fd 2 V + 2M0 t (8^) + - V ' - d r V + Xd r Mg- + Xd, ( -4- ) = xp 2 ^ 



H 2 S 



H 



s 



(17) 



(18) 



where X = FL + M 2 . 

Since the space-time is stationary and invariant by translation in z and rotations in </>, the scalar field can be 
decomposed as 



/oo 
Rkm(t,r) e^ +k ^ dk 
. . - OO 



(19) 



where m is an integer number and k is real. We will be interested in a quasinormal mode decomposition, and therefore 
we also consider a "time-independent" approach, based on a Laplace-like transform 16] 

(20) 

with w extended to the complex plane. Substituting (| 19[) and (|20|) in (p~8|) . we observe that this latter differential 
equation is separable and can be written as 



R{uo,r) = / R km (t,r) e 1 ^ dt 
Jo 



2 H 2 S V H 



X(P 2 



Fm 2 - 2Mmuj - Lui 2 



(21) 
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The evolution of \& in the exterior region of the cosmic cylinders, describing the vacuum, can be calculated straight- 
forwardly. The radial evolution of the massive scalar field in the exterior of a rotating cosmic cylinder is given 
by 



r,d 2 R dR 
x 2 — T +x— + R 

dx z ax 



= 



(22) 



where we defined the coordinate x — (r + ro). Particular cases can be recovered: the static cosmic string and outside 
of the static cosmic cylinder setting ( J = ro = 0) and the rotating cosmic string setting ro = 0. If to 2 > k 2 + (3 2 
equation ([22]) has an exact solution given in terms of Bessel functions. When to 2 < k 2 + (3 2 equation (|22| has an exact 
solution given in terms of the modified Bessel functions decaying exponentially to zero at spatial infinity. Therefore, 
we will analyze only the condition u> 2 > k 2 + /3 2 since, to calculate the quasinormal modes, we want that the field 
behaves as an out-going traveling wave at spatial infinity. 

For the interior of the static cosmic cylinder, the radial evolution of the scalar field has been obtained exactly. 
Applying (|2"Tj) to the respective metric given by ([3]) we have 



d 2 R 



1 dR 



dx 2 tan(:r) dx 



R 



P 2 



sin (a;) 







(23) 



where the radial coordinate x — fo/r has been used. Equation (|23[) has an exact solution in terms of the associated 
Legendre functions. 

For the "flower pot" model the function X and the radial derivative of the determinant of the metric g are 



X = r 2 a\ 



Substituting the functions above in (|2ip the resulting equation reads 



16J 2 



(24) 



, d 2 R 
dr 2 



dR n 
+ r— +R- 

dr 



AJm 



r- 



2 2 

r — to 



= 



with Bessel functions as solutions. 

For the "ballpoint pen" model the function X and the determinant of the metric g are 



X 



im 2 {V\r) 



Using (p?Tj) we obtain 



d 2 R 



dR 



dx 



with the function / given by 



sin {x)— j + sin(x) cos(x)— h R 



(3 2 )sm 2 (x) 



dx 



f = (r — r s ) cos(a;) 



— (m + 2au)fY 



= 



sin(a;) 

~7T 



(25) 



(26) 



(27) 



(28) 



where x = yXr. The equation (|27j) can not be solved exactly. However, an approximate solution has been obtained 
in the quasi-static limit, defined as V\r s <C f- When such a limit is used in (Til)) , the angular momentum and the 
linear mass density become 



J - + 0(x 3 



Xr 2 

^~^±+0(x 3 ) 



(29) 



resulting in an exterior solution that behaves as a quasi-static solution. In this case equation (177)) turns into Bessel 
equation, 



, d 2 R 
dx 2 



dR n 
x— +R 

dx 



(to 2 -k 2 - /3 2 ) 



= 



(30) 
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Another solution of (|27|) can be obtained in the limit of large mass. The necessary condition to be satisfied in this 
limit depends on the interval of the radial coordinate. One can note that in the range < r < r s the function f(r) is 
bounded. Therefore, the limit of large m is reached when 



m max \2aui /Ml 

re[0,r 3 ] 



(31) 



and (1271) becomes 



d 2 R 



dR 



sin (x)— — j + sin(x) cos(x)— hi? 



/3 2 )sin 2 (x) 



(32) 



This equation is solved by associated Legendre functions. 



IV. BOUNDARY CONDITIONS, CAUSALITY, QUASINORMAL MODES 



A. Boundary conditions 



In the static cosmic cylinder backgrounds, the space-time is divided in two different regions. Therefore, we need 
to specify junction and boundary asymptotic conditions. Namely, the scalar field in the interior region \P m must be 
regular at the origin, the field and its derivative must be continuous on the border r s of the cosmic cylinder and only 
outgoing waves can escape to infinity. Hence, these boundary conditions can be summarized as 

\¥ n \ < oo when r -> , (33) 



* m = * > TV = TV when r = r * > ( 34 ) 
ax 11 ox^ 

^out — when r ^ fOQ (35) 

Vr 

where w is a complex frequency of oscillation the scalar field. When the space-time does not have a regular interior 
region, as the static cosmic string, only (|33[) and ([33)) are used as boundary conditions. In the case of rotating space- 
times we need additional conditions to make the quasinormal mode problem well-posed. The difficulty is related to 
the fact that the space-times under consideration are not globally hyperbolic. Although this is not an insuperable 
obstacle, further considerations must be made 

Specifically, there is a region r < r ctc where the temporal coordinate is not globally well defined because the vector 
fields d/dt and d/d(j) are time-like. A compact temporal coordinate (f> leads to closed time-like curves. To deal with 
this problem we follow the approach developed by Novikov, Thorne, Friedman and others, imposing a "self-consistency 
principle" [12| , which roughly states that all events on CTCs influence each other around closed time- like lines in a 
self-adjusted way. One implementation of the self-consistency principle in the dynamical systems discussed in this 
article is done imposing that the propagation of scalar field is periodic in the space-time regions that are identified. 
Labelling these regions by Tj and Tf, we impose "J (r^) = W (T f). In the case of rotating cosmic cylinder space-times, 
this condition is satisfied in the causal region. What we need to do is to extent this condition to the non-causal region. 
We have 

$(t,0,r,z) = V(t,2Tr,r,z) =^ $(0) = $(2tt) . (36) 
With this supplementary condition, the propagation of scalar fields become tractable. 



B. Causality 



In this section an analysis concerning causality is presented. The presence of closed time-like curves can be verified 
analyzing whether < somewhere. In our space-times, the non-causal region are infinitely long hollow cylinders. 

The components of the static cosmic string and cylinder are always positive-definite if fi < 1/4. Thus, the static 
cases can be declared causally well-behaved. The rotating cases need to be analyzed one by one. 
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Rotating cosmic cylinder 





The rotating cosmic string admits a non-causal region with arbitrary size depending on the constants fi and J. In 
Fig. [1] one can see how these constants change the size of non-causal region. There is no restriction for values of J 
but the linear mass density must be bounded. 

The other two cases of rotating cosmic cylinders are more complex than the rotating cosmic string. For simplicity, 
we choose to preserve the interior region, which contains matter, against CTCs, reducing the range of the parameters. 
The exterior of the rotating cosmic cylinder may have CTCs. In case of the "flower pot" model the interior 
component is always positive-definite, and therefore there is no CTC. In the exterior, the condition to have CTCs 
depends essentially on tq being given by 



VTKP 

rctc< (r^)- ro ' with ra = r " =± ^Jf-W- rs ■ (37) 
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If 1 7*0 1 is large enough r ctc < r g resulting in the absence of a non-causal region. In Fig. [2] we show the position of the 
non-causal limit surface r ctc for the rotating cosmic cylinder with the "flower pot" model. 

When the constant 7g" is chosen in (|37|) . the non-causal limit surface always lies before r s , as seen in Fig. [2] (top). In 
this case the exterior of the cosmic cylinder does not have CTCs. On the other hand, when the constant is chosen 
the non-causal limit surface always lies after of r s allowing CTCs. It is shown in Fig. [2] (bottom). This conclusion 
is independent of the value of r s . Therefore, in this case there is no restriction for the value of J. We still have 
(i < 1/4 with the additional constraint \/X r s < n/2. Fig. Q] shows the causal behavior of the interior and exterior 
#00 components for the "flower pot" model. 




FIG. 2: Position of the non-causal limit surface r ctc in the rotating cosmic cylinder with "flower pot" model for (top) and 
Vq (bottom) in (|37[) . The radius of the cosmic string is r s = 5. The dark surface indicates the position of r s and the grey 
surface indicates the position of the non-causal limit surface r c t c . The black arrow shows the position of the non-causal limit 
surface for /i — 1.0 x 10~ 2 and J = 1. 

For the "ballpoint pen" model, both the interior and exterior metrics allow the presence of CTCs. In this model a 
general analysis is quite difficult because of the number of free parameters. Each case must be analyzed independently, 
verifying whether the set of chosen constants satisfies given restrictions. The causal structure this space-time will 
depend on the adjustment of three constants r s , A and a. Analyzing the interior g^ component we can see that it 
can assume negative values. Recalling that we do not want CTCs in the interior solution, the set of constants must 
be chosen such that g^ is positive-definite. Again, the linear mass density must satisfy /i < 1/4 and y/\r s < tt/2. 
The angular momentum is fixed by the choice of r Sl A and a. 

As an example which captures the essence of the argument to be developed here, the set of constants A = 1.0 x I0~ 3 , 
a = 0.2 and r s = 1 keeps the interior g^ > 0, i.e. it is causally well-behaved. Now we need to analyze whether the 
set of chosen parameters allows the presence of CTCs in the exterior. For the exterior of the "ballpoint pen" model, 
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the condition to have CTCs also depends on ro, being given by 



r c tc < 



V16J2 
(1 - 4/i) 



ro 



with ro = r = ± 



y/1 - (1 - 4/g r. 
(1 — 4/i) arccos(l — 4/i) 



(38) 



In Fig. [3] the position of the non-causal limit surface r ctc for the rotating cosmic cylinder with "ballpoint pen" model 
in quasi-static limit is shown. 




FIG. 3: Position of the non-causal limit surface r ctc in the rotating cosmic cylinder with "ballpoint pen" model for r£ (top) 
and Vq (bottom) in (|38[) . The radius of the rotating cosmic cylinder is r 3 = 1. The dark surface indicates the position of 
r s and the grey surface indicates the position of the non-causal limit surface r ctc - The black arrow shows the position of the 
non-causal limit surface for A = 1.0 x 10~ 3 and a — 0.2. 



When the constant is chosen in (|38|) . the position of the non-causal limit surface depends on the value of A. 
For the particular case, A = 10~ 3 and a — 0.2, r ctc is smaller than r s and in the exterior of this cosmic cylinder we 
do not find CTCs. If we choose r = in (|3"8]) the position of the non-causal limit surface also depends of A. But 
in this case r ctc is bigger than r s allowing CTCs. This considerations are general. Constraints on the values of A 
and a results in restrictions for J and /i. Therefore, the angular momentum and linear mass density can not assume 
arbitrary values in the "ballpoint pen" model. 
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C. Quasinormal modes 

In this section we show the exact solutions for the evolution of the scalar field in the cosmic cylinder backgrounds. 
Quasinormal modes are numerically calculated. 

1. Static cosmic cylinder 

The boundary conditions (|3"3"|) and (135[) result in 

R*»(r) = Cx P™(cos(r/f )) , (39) 

{Ci Hi (pr) when uo R > 
(40) 
C 2 H 2 (pr) when w fl < 

where n(n + 1) = (oj 2 - k 2 - (3 2 )f 2 , v 2 = m 2 /(l - 4^) 2 and p 2 = {lj 2 - k 2 - (3 2 ). The functions R in and R out are the 
radial part of the scalar field inside and outside of the cosmic cylinder. One can note that R out depends on the sign 
of the real part of lj, due to the asymptotic behavior of Hankel functions. At the boundary, condition (|34[) implies 



R' m (r s ) _ R' out {n) _ Q 

R m (r s ) R out {r b ) ' 1 ' 

where ' denote a derivative with respect to r. The parameters are related by 

r s =f 6 M , (l-4 M ) = cos(0 M ) , n = r ° sm ^ , p = JL . ( 42 ) 

(1 - 4/x) 87rr<j 

If we choose p and r s , the other constants are fixed and we need only calculating the zeros of (1411) for given k, m, [3 to 
obtain the quasinormal frequencies u>. Some quasinormal modes to static cosmic cylinder were numerically calculated. 
They are shown in Table |TJ 



TABLE I: Quasinormal modes for static cosmic cylinder with the parameters r s — 5.0, p = 1.0 x 10 -3 , fo = 6.3, fi = 0.1, 
rt — 6.4, f3 — 0, k = 0. (Left) Fundamental (n — 1) and high overtone modes (n > 1) for m = 0. (Right) Fundamental mode 
(n — 1) for several values of m. 



n 






UJR + iui 






m 




LOR + iiUI 






1 


-0.66 


- i 


0.54 


+0.66 - 


i 


0.54 





-0.66 - 


i 0.54 


+0.66 - 


i 


0.54 


2 


-1.33 


- i 


0.65 


+1.33 - 


i 


0.65 


1 


-0.21 - 


i 0.42 


+0.21 - 




0.42 


3 


-1.98 


- i 


0.72 


+1.98 - 


i 


0.72 


2 


-0.49 - 


i 0.51 


+0.49 - 


i 


0.51 


4 


-2.62 


- i 


0.78 


+2.62 - 


i 


0.78 


3 


-0.76 - 


i 0.58 


+0.76 - 


i 


0.58 


5 


-3.25 


- i 


0.82 


+3.25 - 


i 


0.82 


4 


-1.01 - 


i 0.63 


+1.01 - 


i 


0.63 


6 


-3.89 


- i 


0.85 


+3.89 - 


i 


0.85 


5 


-1.26 - 


i 0.68 


+1.26 - 




0.68 


7 


-4.52 


- i 


0.88 


+4.52 - 


i 


0.88 


6 


-1.51 - 


i 0.72 


+1.51 - 




0.72 


8 


-5.15 


- i 


0.91 


+5.15 - 


i 


0.91 


7 


-1.76 - 


i 0.76 


+1.76 - 


i 


0.76 



If we display the QNMs in the lji x ujr plane, we can see that the positive and negative lur of the quasinormal 
modes, are distributed symmetrically with respect to the w/-axis. This happens for any value of m. The imaginary 
part, oji, is always negative indicating stability of the static cosmic string against scalar perturbations. The 
oscillates in the four first modes and later grows, making luj more negative. For k ^ and f3 ^ the position of 
modes is pushed away of a^-axis but the quasinormal spectrum qualitative behavior is the same. 



2. Rotating cosmic cylinder 

We now consider the rotating cosmic cylinder with the interior solution of the "flower pot" model. Applying the 
boundary conditions (f3"3"|) , (|34p . (J3SJ) and (|3"o| to the scalar field we have 

R' in (r s ) _ R'° ut (r s ) 

R in (r s ) R°^{r s ) 1 1 
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The functions R m and R out are given by 



R out (r) 



R in {r) = C x J m { Pl {r)) , 

dHl{p{r +r )) when lj r > 
G 2 Hl{p{r + ro )) when < 



(44) 



(45) 



where p 2 = A — 4Jm/r 2 ) 2 — k 2 — /3 2 , m is an integer and ;/ 2 = (m + 4Jw) 2 /(l — 4/x) 2 . Once more the parameters 
associated to the space-time are related by 



16J 2 



ro 



(1 - 4 M )2 



(46) 



If we choose three constants, fi, J and r s , the constant r is fixed and then we need only to calculate the zeros of the 
equation (|43j) for specific values of k,m,/3 to obtain the quasinormal frequencies. 

In Table HT1 some numerically computed frequencies are shown for rotating cosmic cylinder with "flower pot" model. 
In contrast to the static case, here the distribution of uj r is not symmetric with respect to wj-axis. The imaginary 
part u>i can assume negative or positive values depending on the values of in (|46|) . The presence of modes with 
lui > indicates instabilities of the scalar field perturbation. On the other hand, if we choose the constant 7q~ in (|46|) . 
these unstable modes are absent. In this case, the space-time can be considered stable. The relation between and 
the stability issues will be discussed in section IVl 



TABLE II: Quasinormal modes for rotating cosmic string with "flower pot" model with the parameters J = 1.0, /i = 1.0 x 10~ , 
r s = 5.0, /3 = 0, k — 0, ro = —12 (top) and ro = 1.7 (bottom). At the left, fundamental (n = 1) and high overtone modes 
(n > 1) for m — 0. At right, fundamental mode (n = 1) for several values of m. The *'s indicate the unstable modes. 



n 




LUR + i UJI 


m 






LOR+i LUl 




1 


0.064 + 


i 0.041 * 


-0.064 + i 0.041 * 





0.064 + 


i 0.041 * 


-0.064 + 


i 0.041 * 


2 


0.75 - 


i 0.23 


-0.21 + i 0.044 * 


1 


0.026 - 


i 


0.055 


-0.026 -I- 


i 0.034 * 


3 


1.57 - 


i 0.32 


-0.79 + i 0.0082 * 


2 


0.051 - 


- 


i 0.070 * 


-0.56 + 


i 0.026 * 


4 


2.39 - 


i 0.36 


-0.90 + i 0.014 * 


3 


0.15 - 


i 


0.22 


-0.19 + 


i 0.31 * 


5 


3.20 - 


i 0.40 


-1.46 + i 0.031 * 


4 


0.22 - 


i 


0.31 


-0.37 + 


i 0.042 * 


6 


4.01 - 


i 0.43 


-1.57 - i 0.32 


5 


0.29 - 


i 


0.39 


-0.44 + 


i 0.59 * 


n 




ujr + i uji 


m 






LOR + i LUl 




1 


0.46 - ' 


0.45 


-0.46 - i 0.45 





0.46 - 


i 


0.45 


-0.46 - i 


0.45 


2 


1.23 - ' 


0.57 


-1.23 - i 0.57 


1 


0.94 - 


i 


0.39 


-0.83 - i 


0.57 


3 


2.02 - - 


0.63 


-2.02 - i 0.63 


2 


1.43 - 


i 


0.39 


-0.12 - i 


0.15 


4 


2.83 - - 


0.67 


-2.82 - i 0.67 


3 


1.90 - 


i 


0.40 


-0.07 - i 


0.26 


5 


3.63 - ' 


0.70 


-3.63 - i 0.70 


4 


0.30 - 


i 


1.90 


-0.17 - i 


0.30 


6 


4.43 - ' 


0.73 


-4.43 - i 0.73 


5 


0.83 - 


i 


2.09 


-0.26 - i 


0.33 



For the rotating cosmic cylinder with the "ballpoint pen" model, the quasinormal modes have been obtained in two 
different limits, namely quasi-static limit and large m limit. Again the same procedure described before can be used 
to both limits, imposing the same boundary conditions. 

The boundary conditions (f33)) , (f34|) . (f35|) and (f36|) in the quasi-static limit give again the equation (|43[) . The 
functions R m and R out are given by 



R ou \r) 



R m (r) = Ci J m (pi(VXr)) , 

d Hl(p{r + r )) when lo r > 
C 2 Hl(p(r + r )) when io R < 



where now p\ 



/3 2 ) /A. The parameters associated to the space-time are related by 



± VI - (1 - 4/x) 2 r s 

r =r = ±- 



Ar s 2 1 
M~ "IT < 



J-0 



"(l-4^)arccos(l-4^) ' ^ 8^4 

Choosing A and r s , the other constants are fixed and we need only to compute the zeros of (|4"5)l . 



(47) 
(48) 

(49) 
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TABLE III: Quasinormal modes for rotating cosmic cylinder with "ballpoint pen" model in quasi-static limit, with the param- 
eters J = 2.0 x 10~ 5 , £t = 1.0 X 10~ 4 , r s = 1.0, = 0, k = 0, a = 0.2, A = 1.0 x 10~ 3 , r = -2.0 (top) and r = 3.0 x 10~ 4 
(bottom). At left, fundamental (n = 1) and high overtone modes (n > 1) for m = 0. At right, fundamental mode (n = 1) for 
several values of m. The *'s indicate the unstable modes. 



n 




WiJ + i LOI 




m 


LOR - 


-iuj 




1 


0.30 + i 0.41 


-0.30 + 


i 0.41 * 





0.30 + i 0.41 


-0.30 + 


i 0.41 * 


2 


2.96 - i 0.94 


-2.96 - i 


0.94 


1 


1.84 + Z 0.25 


-1.84 + 


i 0.25 * 


3 


6.16 - i 1.28 


-3.83 + 


i 0.18 * 


2 


0.34 + i 0.47 * 


-3.06 + 


i 0.27 * 


4 


9.33 - i 1.47 


-5.52 + 


i 0.17 * 


3 


1.10 + i 0.70 * 


-1.11 + 


i 0.70 * 


5 


12.5 - i 1.62 


-6.16 - i 


1.28 


4 


1.92 + i 0.86 * 


-1.92 + 


i 0.86 * 


6 


15.6 - i 1.73 


-7.01 + 


i 0.17 * 


5 


2.76 + i 0.99 * 


-2.76 + 


i 0.98 * 


n 




LOR + i LOI 




m 


LOR - 


- i uj 




1 


2.61 - i 5.54 


-2.61 - i 


5.54 





2.61 - i 5.54 


-2.61 - 


i 5.54 


2 


5.94 - i 5.72 


-5.94 - i 


5.72 


1 


1.07 - i 5.47 


-0.66 - 


i 5.46 


3 


9.17 - z 5.87 


-9.17 - i 


5.87 


2 


2.74 - i 5.54 


-2.21 - 


i 5.60 


4 


12.4 - i 6.00 


-12.4 - i 


6.00 


3 


1.00 - i 5.86 


-0.53 - 


i 5.90 


5 


15.5 - i 6.10 


-15.5 - i 


6.10 


4 


2.37 - i 6.15 


-1.94 - 


i 6.27 


6 


18.7 - i 6.18 


-18.7 - i 


6.18 


5 


0.85 - i 6.77 


-0.49 - 


t 6.81 



TABLE IV: Quasinormal modes for rotating cosmic cylinder with "ballpoint pen" model in large m limit with the parameters 
J = 0.1, n = 1.0 X 10" 2 , r 3 = 10, (3 = 0, fe = 0, a = 1.0, A = 1.0 X 10" 3 , r = -20 (top) and r = 0.4 (bottom). At left, 
fundamental (n = 1) and high overtone modes (n > 1) for m = 70. At right, fundamental mode (n = 1) for several values of 
m. The *'s indicate the unstable modes and f's indicate the purely real modes. 



n 




LOR +i LOI 


m 






1 


2.20 T 




-2.20 T 


60 


1.88 T 


-1.88 T 


2 


2.47 + i 3.99 




-2.31 + i 4.30 * 


62 


1.94 f 


-1.94 f 


3 


2.80 + i 3.86 




-2.84 + i 4.08 * 


65 


2.04 f 


-2.04 f 


4 


2.96 + i 3.78 




-3.02 + i 3.99 * 


67 


2.10 f 


-2.10 f 


5 


3.12 + i 3.70 


* 


-3.20 + i 3.90 * 


70 


2.20 f 


-2.20 f 


6 


3.29 + i 3.61 




-3.38 + i 3.80 * 


72 


2.26 f 


-2.26 t 


n 




LOR +1 LOI 


m 




LOR + i LOI 


1 


0.033 - i 5.00 




-0.15 - i 4.70 


60 


0.18 - i 4.30 


-0.10 - i 4.05 


2 


0.21 - i 5.00 




-0.32 - z 4.70 


62 


0.19 - i 4.44 


-0.11 - i 4.18 


3 


0.39 - i 4.99 




-0.48 - i 4.69 


65 


0.11 - i 4.65 


-0.045 - i 4.38 


4 


0.57 - i 4.98 




-0.64 - i 4.67 


67 


0.11 - i 4.79 


-0.057 - i 4.51 


5 


0.75 - i 4.96 




-0.81 - i 4.66 


70 


0.033 - i 5.00 


-0.16 - i 4.71 


6 


0.93 - i 4.94 




-0.97 - i 4.63 


72 


0.039 - i 5.14 


-0.17 - i 4.84 



Quasinormal modes for the "ballpoint pen" model in the quasi-static limit were obtained numerically being shown 
in the Table Mil The imaginary part u>i can assume negative or positive values depending on the values of r in 
(|4"5)) . as in the "flower pot" model. In the stable case, the first modes have higher loj than u>r. This indicates that 
these modes are strongly damped. In the same way the presence of u>i > would indicate instabilities of the scalar 
perturbation. However, if we choose Tq in (|49|) these modes are absent. 

Using the boundary conditions (f3"3"l) , ([M]) . and ([5r])) in the large m limit, equation (HS")) is obtained with 

i? m (r) = Ci P^(cos(V^r)) , (50) 

CtHlipir +r )) when > 
i?°"'(r) = <; , (51) 

C 2 Hl(p(r + r )) when w fl < 

where n(n+ 1) = (cj 2 — fc 2 — /3 2 ) /A. In Table [TVl some quasinormal modes for the rotating cosmic cylinder "ballpoint 
pen" model in large m limit are listed. The sign of the imaginary part uij depends again on the values of r^. The 
presence of loj > indicates instabilities of the scalar perturbation. On the other hand, if we choose in ([4^)) these 
modes are absent and we can declare the stability of the space-time under scalar perturbations. The precision of the 
approximation is higher in the first overtones. In this limit the field presents purely real modes of oscillation but they 
arise only when instabilities are also observed. 
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3. Static and rotating cosmic string 

In the case of static and rotating string, the application of the boundary conditions (|34| and (|35[) lead to Hankel 
functions. However when regularity condition at the origin is imposed, both Hankel functions diverge. Hence these 
space-times do not display quasinormal modes of oscillation. 

V. FINAL REMARKS 

The results in this paper concern the stability of cosmic cylinder space-times and its relation with closed time-like 
curves. The data obtained in previous sections show that, when the parameters are chosen in a such way that the 
scalar quasinormal modes have always uj < then an inspection of the causal behavior reveals the absence of CTCs. 
On the other hand, if the parameters are chosen to guarantee the presence of CTCs in the space-time then the scalar 
quasinormal modes will present unstable modes u>i > 0. 

For the static cases and for the rotating cosmic strings and cylinders it is observed that the solutions can be divided 
in two groups: the "group S" , where Tq is chosen in Eqs. ([38 |) . (l46p . (|49|) : and the "group I" , where Tq is chosen. Thus 
the constant Tq controls the appearance of CTCs in these cosmic cylinder space-times. The results can be stated 
schematically as 



Scalar stability 
toi < 


Existence of CTCs 
94><t> < 






group S: 
r+ in Eqs. (JMD , (© , © 


group I: 
ro in Eqs. (US]) , flU]) , @g) 






Absence of CTCs 

9<t><t> > 


Scalar instability 

LUI > 



Summarizing our results we conclude that cosmic cylinder space-times that present closed time-like curves are 
unstable against scalar perturbations. Since the scalar field propagation gives us strong signals of instabilities in the 
considered space-times when the CTCs are present, our results for the scalar perturbation lead us to conjecture that 
the presence of closed time-like curves is one source of space-time instabilities. The demonstration of this conjecture 
could open a new window in the investigation of causality. 
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